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Scalar-tensor propagation of light in the inner solar system at the millimetric level 
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In a recent paper [l[, motivated by forthcoming space experiments involving propagation of light 
in the Solar System, we have proposed an extention of the IAU metric equations at the c -4 level in 
General Relativity. However, scalar-tensor theories may induce corrections numerically comparable 
to the c -4 general relativistic terms. Accordingly, one first proposes in this paper an extension 
of pi to the scalar-tensor case. The case of a hierarchized system (such as the Solar system) is 
emphasized. In this case, the relevant metric solution is proposed. Then, the geodesic solution 
relevant for propagation of light in the inner solar system at the millimetric level is given in explicit 
form. 
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I. INTRODUCTION 



Forthcoming space missions and missions in project - such as LATOR [Hj], TIPO [13 , ASTROD [T|], PLR [Ig|, 
ODYSSEY [201 or SAGAS [2l| - will require distance measurements at millimetric level in the Solar System. This 
corresponds to time transfer at the precision 10~ n s. As argued in |l|, this requires a complete Solar System metric 
at the c -4 level, in order to describe the laser links involved in such experiments. This has been proposed in the 
framework of the General Relativity (GR) theory in [l|, leading to an appropriate extention of the metric equations 
| ■ recommanded by the IAU2000 resolution [2| . 

C_) | The relative amplitude of the relativistic effects is of the order of e = GM/rc 2 , where M and r are some characteristic 

<sj i mass and distance. The so-called first order terms are of order e, the second order terms of order e 2 . In the inner 
5_i ' Solar System, e is typically of the order 10 -8 , and can be sensitively greater (10 -7 ) for photons entering well inside 
&Q, Mercury's orbit, and can even be as large as 10~ 6 for photons grazing the Sun. 

On the other hand, there is a surge of interest in scalar-tensor (ST) theories since about two decades. Indeed, 

I , the gravitational sector of a lot of tentative fundamental theories, like string or (modern) Kaluza-Klein, turns to be 

^. ■ described by a metric tensor plus a scalar field (Brans-Dicke [3, [I| or not Q). Besides, alternative theories to GR 

C\) | (including ST) are also sometimes required by some authors to deal with the so-called dark energy (cosmological 

r ■ level), dark matter (galactic level), Pioneer or fly-by anomalies (solar system level) problems. Phenomcnologically, 

the divergence between ST theories and GR is quantified by the Post-Newtonian (PN) 7 factor (= 1 in GR) entering 

the c -2 term in the space-space components of the metric tensor. From the present observations, I7 — 1| can at best 

C ' reach values of the order of 10 -5 [7[. But, some theoretical considerations strongly suggest 7 could have been driven 

from any " initial" value to a value close to unity by the cosmological expansion (more precisely, ST theories are driven 

to GR, as soon as these theories fulfil some (not very constraining) conditions), and even claim the 1 — 7 current 

value should be of order 10~ 7 or 10 -8 [3 19|. On the other hand, c -4 terms are typically 10 6 to 10 8 smaller than c~ 2 

terms. All these reasons make relevant the extension of the gravitational framework proposed in [l[ to encompass the 

ST case, as much as c space-space metric terms have to be taken into account in light propagation problems. 

In principle, this would require new definitions of multipolar moments at the c -2 level. But other publications 
proposed c -2 multipolar moments in ST theories [10j | or in a parametrized post-newtonian framework [11]. Thus we 
do not discuss this point in this paper and focus only on the c -4 metric side of the problem. Another reason to 
discard this point is that non-monopolar terms turn out to have numerically negligible contributions for links we are 
interested in in this paper (inner Solar system case). 

In section [Hi one defines a terminology relevant to the considered problem. Since the Einstein conformal representa- 
tion plays a central role in the approach followed in this paper, the section UlTl is devoted to the conformal link between 
the representations of ST theories, and to the related notations we will use. Section HVl is dedicated to the derivation 
of the ST field equations up to O (c~ 5 ) terms. In section[V] these field equations are rewritten up to O (c~ 5 , cj _1 c^ 4 ) 
terms, for applications taking explicitely our present experimental knowledge on gravity into account. Considering 
applications to Solar System-like systems, one defines hierarchized systems in section [VT] This is made by defining a 
quantity \x that quantifies how much the system is gravitationally dominated by its most massive body. Accordingly, 
the field equations are rewritten up to O (c~ 5 ,u!q c -4 , ^c~ 4 ) terms. In this case, the explicit metric solution that can 
be used in relevant applications is written in well-suited coordinates. Finally, in section VII, we give the explicit form 
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of isotropic geodesies relevant for time transfer and ranging problems in the inner solar system at the millimetric level 
of accuracy (required by forthcoming space missions or missions in project). 

II. TERMINOLOGY : DEFINITION OF THE PN/BM AND PN/RM METRICS 

The PN approximation is based on the assumption of a weak gravitational held and low velocities for both the 
sources and the (test) body (ie. velocities of the order yjGM/r or less, M being some caracteristic mass of the system). 
It formally consists in looking for solutions under the form of an expansion in powers of 1/c. The usually so-called 
nPN order terms in the metric, leading to c~ 2n terms in the equation of motion of a body describing a bounded orbit, 
are terms of orders c~ 2 ™~ 2 in goo, c~ 2n ~ 1 in goi and c~ 2n in gij. In this case, the Ricci tensor components have to 
be developped the same way as (18) in |l|. In the present paper, a metric developped this way will be refered as the 
nPN/BM metric (BM meaning "Bounded Motion" for test particles). It is particularly well-adapted for studying 
bounded motions in systems made by non-relativistic massive bodies, as the Solar System is. 

However, since we are interested in the propagation of light, we are led to relax the hypothesis on the velocity of the 
test particle whose motion is considered. Of course, this doesn't change the metric, but the terms to be considered in 
the metric components are not the same as in the PN/BM problem. Indeed, the terms leading to c~ 2n terms in the 
equation of motion of a test particle moving at relativistic velocity are terms of order c~ 2n in both goo, g^i and gij. 
In this case, the Ricci tensor components have to be developped the same way as (19) in [if. In the present paper, 
a metric developped this way will be refered as the nPN/RM metric (RM meaning "Relativistic Motion" for test 
particles). A PN/RM metric is particularly well-suited for studying relativistic motions of test bodies (for instance, 
the propagation of light) in systems made by non-relativistic bodies, as the Solar System is. 

The present paper deals with the PN/RM problem since we are concerned in propagation of light. 



III. THE ST THEORIES IN EINSTEIN VS JORDAN REPRESENTATIONS 



The Jordan representation of the ST theories is described by the action 
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In this representation, the gravitational sector of the theory is described by the Jordan metric g a g and the scalar 
field $, while the non- gravitational fields are symbolically represented by *&. The scalar field couples directly with the 
metric, leading to rather complicated field equations. Besides, the kinetic term associated to the scalar field doesn't 
have the standard form, involving a scalar field function u>, the function characterizing the ST theory we are dealing 
with. On the other hand, the non-gravitational lagrangian Lng doesn't depend on the scalar field, leading to simple 
equations of motion (V„T a ' 5 = 0), with the nice consequence that the weak equivalence principle applies in this 
representation of the theory. 

To overcome the just mentionned drawbacks, one could be tempted to resort to a dependent variables change 
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chosen in such a way that the action ([lj transforms into 
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Lng depending on \l/, g and ip in a way to be precised later (R and g correspond to R and g, but with g a p replaced 
by g a a). Since the scalar field doesn't couple with the metric, ([2]) is refered to as the Einstein representation of the 
theory. 

The form ((2|) is achieved by considering a conformal transformation of the metric 



g a/ 3 =A(tp) 2 g a p. 
From the induced transformation of the Ricci scalar [12| , and up to a divergence term, $Z§ turns into 
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Comparing Q with ([T]) suggests : 

- the link between the Jordan and Einstein representations of the scalar field 



* = ?; ( 5 ) 
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the Einstein representation of the non gravitational lagrangian 

L NG = A(ip) 4 L N G{^,g,iu) 
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(HJ) and (|6]) lead to the link between the functions w ($) and A (<p) (equivalently characterizing the considered ST 
theory) 

([7]) requires w > —3/2. This results from the fact one has imposed the sign of the scalar kinetic energy term in Einstein 
representation (|5J| in order to ensure the dynamical stability of the theory [13J . 

The link between the stress tensor components in the two representations follows from the general stress tensor 
definition 



5g a P — -> 6 / d^Xy^jL 
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Sg a — > S / d 4 x\/^L NG = — - / d i x\/^f a p6g a . 

Since \j —gL^G = \J '—gL^c and since the scalar field is not varied in this metric variation process (no ambiguity 
since the two versions tp and $ of the scalar field are related in the non metric dependent way (J5j)) , it directly turns 
that 

T a p = A T a p. 

For the mixed and contravariant components, it follows 

T^ = A*t£ (=>f = A i T), f a/3 = A 6 T a P 
the indexes being raised/lowered by the metric involved in the corresponding representation. 

Eliminating A 

It is clear A can be eliminated between the two representations of the scalar field using (JSJ . This way, the considered 
ST theory is represented by the function 4> (tp) in Einstein representation. From ([7]), this function is linked to u ($) 
by 

rfln$ x 2 



^ + 3) {—) = 4 - (8) 

The conformal transformation (j3]), the link between the non gravitational lagrangians and stress tensors now write 

g a p = $5a/3 (9) 

L NG = 1>- 2 L NG (#, $ _1 <7 M „) 



T a /3 = $ T a p, 

T = $- 3 T a/3 . (10) 

As it turns from equation (|8]l. tp is denned up to a sign and an additive constant. When needed in the following, 
the sign will be fixed by the choice 

V2^T3^ = 2. (11) 

dp 

IV. THE FIELD EQUATIONS 

A. Using the Einstein representation gravitational field variables 

From ©, the field equations can be written 

n I i — -7\ a \ 27rG- / — rdln$ . „. 

the function $ (p) characterizing the ST theory explicitly entering the scalar field equation. As usual in PN approxi- 
mation, let us write the scalar field as 

(2) (4) 

<P = ^0 + ^ + ^ + O(c- 5 ) (13) 

(2) (4) 

where p$ is constant and p and p are zeroth order terms. (Remark that, since p is defined up to an additive 
constant, it is not restrictive to set p$ = 0.) As a consequence, it turns out that, under the standard PN assumptions, 

R = O (c -4 ) , so that the Strong Spatial Isotropy Condition (SSIC) [1J] applies in this representation. It is then 
possible to choose a coordinate system in which the Einstein metric takes the following form at the 2PN/RM level [1| 

.9oo = -l + ^-^ + O^ ( 14 ) 
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Putting &((fo) — 1 is not restrictive since $ enters (fT2")l throught its logarithm derivative. Hence, the scalar field 
function develops as 

$ (p) = 1 + ^*$ + 1 (*$ + ^f/^j + O (c- 5 ) (15) 

where $q and $q stand for the values of the derivatives of $ at po- Setting 



]_ ^ T oo + T kk^ ( 16 ) 
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a* = -T 0i 
c 

a ij = rpij _ T kk § .. ( == ^ a kk = _ 2T kk ) 

(which, from standard PN assumptions, are c order quantities) and using (fT0|). the (00), (Oi) and («j) field equations 
(|T2"j) lead respectively to 

Auj + -i (3a«tfl + Ad tk w k ) + ^<P' { pAw = -4^Gct + O (c" 3 ) (17) 



5 
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®ij (m) = diwdjw - d t (d t Wj + djWi) - 2S ij d t {d t w + d k w k ) + ^■kGcj 13 +d l ipd j ip +0 (c" 1 ) (19) 

where &ij is defined, as in [lj, by 

&ij (tm) = d ik Tj k + d jk r ik — Ary — dijT kk . 
The scalar field equation gives 






A$ + ^ -d t $ + A^\+^(W -^-\ #A? = -2^ a + ^- ) + O (c^) . (20) 
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Remark that, in contrast to the GR case, the (00) equation is not linear, because of the c -2 ip Aw term. The scalar 
field equation also contains a non-linear c -2 term. 
Now, combining (|T7|) and (f2"0")) leads to 



Accordingly, let us choose 



a( ( J-1$>) =o(c-i). 



% = U>' Q w. (21) 



(4) 

Hence, defining x = ¥> /$o 

^ = ^o + ||^ + §x + 0(c- 5 ). (22) 

The metric field variables iu, u>^ and Ty are now decoupled from the scalar field \. The system constraining w, W{ 
and nj now writes 

Aw + 1 (35 ttW + 4d tfe w fe ) + ^ fcfwAw = -47rGa + O (cr 3 ) (23) 

Aw, - d lk w k - d ti w = -AirGa 1 + O (c~ 2 ) (24) 

9y (r«) = -d 4i u>j - ftjWj + ( 1 + J$o 2 J ftt«S|,-t(; - 2dij (d tt w + d tk w k ) + AnGa^ + O (cr 1 ) (25) 

X being obtained in a second step, by solving 

A X - 2 {d tt w + d tk w k ) + I « - $£) wAw = -2wGa kk + O (c" 1 ) . (26) 

Remark that Aw may be replaced by — AirGa in the non linear terms of equations (|23j) and (|20l) . 

B. Back to Jordan representation 

Let us use the function w ($) and its derivative u/ ($) instead of 3>' (<p) and $" (<^). One finds, using (ITT1) 

o 
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One now goes back to Jordan representation using (|9]), with, from (fT5|) and (|22l) . 
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Now let us put 
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and let us define 



(U,Ui, Uij , P) = — (w,w it nj , x) 



and the related quantities (W, W», Wy) by 
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Using (fT4")l and (J2"T)) . one gets the Jordan metric 
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where the functions (W, Wj, Wjj,P) satisfy the following field equations - after some algebra from (|23H26p and (|16p 
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In (|2"9")l . one has set 



and, for the matter part of the equations 



J = d t U + d k U k 

= d t w + d k w k + ( c - 2 ) 
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Let us remark that the quantity 2(/? — 1)/ (1 — 7) (= uj' q (2ojo + 3) (2ujq + 4)~ ) is not diverging, even if 7 is 
(arbitrarily) close to unity. Besides, no new PN parameter appears neither in the c space-space part of the metric 
nor in the corresponding field equations, as stressed in [6(. 

This form is relevant in all sufficiently weak gravitational field, even in systems where the ST theory is not very 
close to GR, i.e. where the PN parameters 7 and /3 are not close to unity. A priori, this may occur even if 7 and /3 
are close to unity in some (other) regions of the universe, as in the Solar System, as soon as the ST theory is not (in 
some sense) close to the Brans-Dickc one (in Brans-Dicke gravity, to doesn't depend on the scalar field, so that it has 
the same value in all the space-time regions of the universe). 

Let us point out that the numerical values of the coefficients j 2 + f3 — 1 and 7 + 1 entering g^ in (29) can be chosen 
independently one to another, since both j3 and 7 enter these coefficients. Hence, there is no a priori relation between 
the coefficients of the c -terms W 2 and W^ , contrary to what may be suggested by the form of the metric chosen 
in [34j (in the one-mass case). More precisely, this doesn't mean the form chosen by [34[ is uncorrect, but rather that 
this form doesn't encompass the (general) ST case (but it encompasses the GR case, as it must be). 

C. Harmonic gauges 

Since the use of the harmonic gauge (HG) is recommended by the IAU, let us consider the field equations in this 
gauge. Of course one has to specify the representation in which the HG is prescribed. The Jordan HG condition reads 

9 L a/3 — U 

and it leads to, for the space (a = k) component 

{ 1 -l)d k U = 0(c~ 2 ). (32) 

As expected from known results in GR J14j . this condition reduces to a triviality in the case 7 = 1. On the other hand, 
if 7 7^ 1, (|32t shows that the coordinate system in which the metric takes the (Jordan) form (f28|) , corresponding to 
SSIC in Einstein representation, doesn't encompass (Jordan) harmonic coordinates in the ST case. In other terms, 
(Jordan) harmonic coordinates are incompatible with the SSIC in Einstein representation. 
One could rather choose to impose the HG condition on the metric in Einstein representation 

9 r a/J = 



since the Einstein metric (|14|) satisfies the SSIC. From (|T4j) . this means one imposes w, Wi and 77, to satisfy 

d t w + d k w k = 0(c~ 2 ) 
d k T lk -\ dl T kk+ d tWl = 0(c-i). 

Translated in terms of (U, Ui, £%), this takes exactly the same form, i.e., using (|30II31[) 

J = 0(c- 2 ) , J i = 0(c~ 1 ). (33) 

It turns out this corresponds to the Nutku gauge constraints [10|, LL5( , meaning that imposing the HG in the Einstein 
representation is equivalent to impose the Nutku gauge in the Jordan representation. Using (|33)) . the three first 
equations of (|2"9")l take the reduced form 

BW + 1 + 2 ^^WAW 

= -47rG e// £ + O (cr 3 ) (34) 

AWi = -AirGeffZ 1 + O (c~ 2 ) (35) 

AWij + diWdjW + 2(1-/3) 6 tJ WAW = -47rG e// £ lJ + O (c^ 1 ) (36) 
while the fourth equation of (|2"9"|) and the harmonic constraints read, using (|30p and (pTTj) 
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AP + 2- WAW = -4irG eff - + O (cr 1 ) (37) 

d t W + d k W k = O (c~ 2 ) (38) 

d k W lk - \d,W kk + d t Wi - ^-^P = O (cr 1 ) (39) 



or equivalcntly (after elimination of the scalar field P) 

d t W + d k W k = O (c~ 2 ) (40) 

d lk W lk - -AW kk + duWi + 09 - 1) WAW = -2nG e ff^— ^S fcfe + O (c" 1 ) (41) 

OtkWjk + BtiWj = d jk W ik + d t3 W t + O (c- 1 ) . (42) 

The last equation refers to the fact that diP (given by ([Ml) - ) is a gradient. 

Note these equations are coherent with 1.5PN/BM equations assumed in [llj . 

V. RELEVANT FIELD EQUATIONS CONSIDERING PRESENT CONSTRAINTS ON GRAVITY 

A. Without making the HG choice 

In the inner solar system, the gravitational field is such that 

% ~ 10- 6 to 1(T 8 . 

On the other hand, from experimental/observational constraints [7( 

[-y — 1| < 1(T 5 ic. uj > 10 5 . 

This means 7 — 1 (or u>q ) could be considered numerically as a c _1 (at best) order quantity. Hence, it is convenient 
to present the metric under the form of a generalized development in both powers of c~ 1 and loq . The useful metric 
resulting from (|2"5)l reads (if u>' is not "unreasonably large") 
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where W, W* and W,j satisfy, from (|29l) 

DM/ + -d t J = -A-KG eff o + O (c-V^c- 2 ) 



AWi - <% J = -47rG e// cr i + O (c~ 2 ) 

AW W + diWdjW - d.J, - djJi - 25 tJ d t J = -4TTG eff a^ + O (c' 1 ,^ 1 ) 

and where Ji reduces to 

Ji = d k W lk - ^diW kk + d t Wi + O (uj^ 1 ) . (43) 

One remarks the field equations take exactly the same form as the GR case [l| (with G replaced by G e //). The only 
remaining reference to the scalar field is reduced to the presence of the 7 PN coefficient in the metric tensor. Related 
to this, the field equation on P is dropped out. 

B. Making the HG choice 

The corresponding harmonic equations to be used when considering known constraints on loq reads 

DW = -4irG eff a + O (c~ 3 , w^V 2 ) 

AWi = -47rG e// cr 4 + O (c- 2 ) (44) 

AWij + diWdjW = ~A-KG eff a lj + O (c~\ c^ 1 ) 



with gauge conditions 



d t W + d k W k = O (c- 2 ) (45) 

d k W lk - \d % W kk + dtW, = O (c" 1 ,^ 1 ) . (46) 



VI. EXPLICIT HARMONIC METRIC RELEVANT FOR HIERARCHIZED SYSTEMS 

A. Hierarchized systems 

Let us consider the case where the system is composed by bodies of masses Ma- Let us consider one of these bodies, 
named S, of mass Mg. Let us define the parameter 

One defines a hierarchized system as a system in which the body S can be chosen in such a way that 

H <C 1. 

In such a system, the body S will be hereafter referred as the " star" , while the other bodies will be referred as the 
"planets". 

In the general relativistic iV-body problem, multipolar moments of a body A are defined in the coordinate system in 
which this body is, in some sense, at rest. These moments arc affected by coordinate transforms through a "Lorcntz- 
like length contraction effect" . These effects being of order (u/c) , where u is the relative velocity between the two 
frames, the induced effects in the metric components are of order c -4 , since potentials are at least c -2 terms. 

In hierarchized systems, the velocity of the body S is of the order of 



GM S 

VS ~ fJ-VB ~ fi\ 

V r S-B 

where B is the most massive planet (and r$-B the distance between B and the star). All the Lorentz-like contraction 
terms have a form like 

GM A v\ 
re 2 c 2 

If A is a planet (A ^ S), this term is at best of order O (^c -4 ), since Ma < nMs- If A is the star (A = 5), this term 
is of order O (vgC 4 ), ie. O (/i 2 c -4 ). Hence all these terms are, at least, of order O (/ic~ 4 ). 

Let us also point out that, since at this level the metric depends on time through the positions of the star and the 
planets only, all the terms containing the operator dt are at least of order 0(/j,). Hence, equations (|4"4"]) -([4l) ]) lead to 

AW = -AwGeffa + O (c- 3 ,coo l c~ 2 , nc~ 2 ) 

AWi = -4TrG e ffa l + O (cr 2 ,^- 1 ) (47) 

AWij + dtWdjW = ~4irG eff a l: > + O (c -1 ,^ -1 ,/*) 

with gauge conditions 

d t W + d k W k = 0(c- 2 ,nc- 1 ) (48) 

d k W ik - ^diW kk = O (c-\ Uo\n) . (49) 

Related to this, A^ defined in [2( leads to numerically negligible terms (see (11.4.8) in [HJ for the ST version). 
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B. Application to the solar system 

In the Solar system, the most massive body S is the Sun and one has 

H~ 10~ 3 . 

Thus, it is legitimate to consider the Solar system as a hierarchized system. Note that, at best, only J 2 , J4 and Jq 
planetary multipolar coefficients (for giant planets) could have a significant impact on laser ranging experiments at 
the required accuracy (see J22M25I ] for giant planets' multipole moments values). Hence, taking advantage that the 
Solar multipolar terms are very weak, the solution of the field equations (|47j). suitable for millimetric accuracy in 
propagation of light, with the harmonic constraints (|48H49[) given in barycentric coordinates turns to be 

900 = -l + ^[W (i,x) + Wi(t,x)]--^ + O(c- 5 ,a; - 1 c- 4 ,Mc- 4 ! J 2 s c- 4 ) (50) 

SW = -2^±±W\t 1 x)+0{c-\^\jlc- A ) (51) 



c° 



/ 9 2W^ \ W^ 

9ij =1 + 4 [Wo(t, x) + W L (t, x)} + — f° % +4^ + (c- 5 ,^ 1 ^ 4 , [ic-\ J 2 s c- 4 ) (52) 



where 



W (t,x) = y W Afi , and W A ,o = G eff ™ A . (53) 

3 R 2n ( k a ■ r a\ 

W L {t,x) = yW A , L , with Wa,l = -G eff J2 m aJL ^TT p 2„ -*—*- ■ (54) 

a 1 T a \ 7 .4/ 



W l {t, x) = J2 W\(t, x), with W A (t, x) = G e 



If 



2)i 



^a(i + ^ rAR 2 AD 

r A \ 

1 (n . M sY ((■ -4)(^-4) 



\ l + T. J 2nJn^ 



wm - -A G *ff rsJ I 




where one has put 



r A (t,x) = i — 4^) and rA(*,a?) = |?^4(i,x)| = y (x l - x A )(x l - x\). 



M A , r a , v A and S A are repectively the mass, the position and the velocity in barycentric coordinates, and the total 
angular momentum of the body A. R A and J„ are the radius and the mass multipole coefficients of the body A. P n 
are the Legendre polynomials and k A denotes the unit vector along the local Z A axis of each body A. The differences 
with the IAU2000 metric [2[ lie in the presence of both the PN parameter 7 and the c -4 space-space metric term. 
The multipolar term Wl in gij that has been neglected in the IAU2000 metric - thanks to numerical considerations 
in the 1PN/BM case - has to be considered here as well. Accordingly, multipolar terms enter also the time-space 
component of the metric (<7oi)> an< ^ could lead to measurable effects, depending on J2 n orders of magnitude. 

While the rotational term in the time-space component of the metric is given as the usual Lense-Thirring term, 
slight modifications (spin multipoles) can in principle appear due to the differential rotation of the bodies. However, 
Solar seismology suggests (26| that the Sun's tachoclyne is at about 0.7 Sun radius. Then, the mass concerned by 
the differential rotation is of order of a few percent of the total mass and thus, it might not lead to measurable 
effects. But, incidentally, time transfer and laser ranging experiments could suggest a new way (independent of Solar 
seismology) to test our knowledges on the solar interior dynamics. Another point coming from the solar seismology is 
that the solar core (r < 0.2 Sun radius) may rotate faster than the external layers [28l - l3(j |. Since, the core represents 
a great amount of the total mass, this could affect the propagation of light at a level depending on the total angular 
momentum value (<Ss). Depending on the solar internal structure model, this may happen at the millimetric level. 
Using a simplified model, we show in appendix lAl how the non-rigidity of the rotation affects the metric. 
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Let us note that the Wl terms can be neglected for inner solar system millimctric laser ranging experiments, such 
as Mars laser ranging for instance. 

We emphasize again that neiher the /3 parameter nor any e parameter are required (e - corresponding, sometimes 
to A [3l| or 5 [lj| [33[ - being some PN parameter often considered in the c -4 space-space metric term 32|, [34| ) . This 
is because both the former and the latter give too small deviations from GR to be considered in c -4 Solar system 
photon's trajectories calculations. This fact is known for the former from [7j and is then obvious for the latter since 
it is a function of 7 and /3 in the (non-massive-) ST theories considered here (as expressed in equation (J28D ). 



VII. ISOTROPIC GEODESIC SOLUTIONS RELEVANT IN THE INNER SOLAR SYSTEM 

A. The metric to be used for time transfer at the millimetric level 

Space missions like LATOR [H[, TIPO [l7|, ASTROD [p, Phobos Laser Ranging [3 require laser links at the 
millimetric level in the inner solar system. As discussed in PI , a full c -4 metric like (|50M52|) is needed, and the transfer 
requires writting the isotropic geodesic equations attached to this metric. However, it turns that some terms in the 
list (|53M56[) can be neglected. Indeed, an order of magnitude analysis (relevant for propagation of light restricted to 
the inner solar system) shows the following terms are not of numerical relevance : 

• non-monopolar terms {Wl and j£ n terms in W 1 ) ; 

• planetary spin terms (S terms in Wp lanets ). 

Hence, the metric components we consider are restricted to (with Sg = (constant) solar moments) 



.9oo = 


1 (2) (4) 

= -1+ .9 00+ .9 00 


got = 


(3) 

= 9oi 


9ij = 


= °ij + .9 ij + 9 ij 



(57) 



with 



(2) 

9 00 



(4) 

9 00 

(3) 

9 ot 



(2) 

9ij 

(4) 

9ij 



t-^2GM A 1 _ W2) 
Z^ a — -1^9 (A)00 



(2 

5 00 



i\\ 



1 2 



n G _i -De ,. , v— \ n ; 2GMa 1 

-2(l + 7)-3^5|-f-(l + 7 )E^^^- 
c r s A c ta 



x (2) 
10ij .9 00 

A (4) ^C 
~°ij 9 00 + Sij 



with e« = — T^ ^T^ - 3 



where D l A (t, x l ) = x l — x l A (t) (components of ta {t, x)), r A = D l A D l A and /3 A = v A /c. 

(4) (4) 

The calculations have been made with a tracking coefficient e in front of g 00 and g ,_■ in order to separate terms 
comming from the fourth order part of the metric from those comming from (2nd order)- (2nd order) geodesic equation 
coupling terms effects in the geodesic solution. However, for convenience, one has put it to unity in the following. 



One gets 
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oo 



(3) 
pO _ pO 
1 oo — l 00 



(2) 
"pA: 

— l 00 

(2) 

— l Oi 

(3) 

pfc 

^ L Oi 

(3) 

= r° 






•J 



(2) 



(4) 
p/c 
1 00 

(4) 

r° 



(4) 
y-ifc 



with, taking into account only relevant terms 



(2) 



Oi 



(2) 
pfc 

J-oo 

(2) 



(3) 

-Lqo 



(3) 

it 



(3) 
p/c 
1 Oi 

(4) 
1 Oi 

(4) 

"pfc 

i 00 

(4) 

"p K 



(2) 

= Zi 

(2) 

= z k 



(2) 



(2) 



(2) 



= -jdik Zj- ySjk Z j + ydij Z k 

v- / (2) 



(A) 



(3) 

Zi 



(3) 



(A)l 

{A) 



(2) 



(A)l 



(3) (3) 

= — Z ik + Z ki 



l&ik^2 PAe Z {A)l 
(A) 



(2) (2) 



■2g 

(4) 

- z ijk 



(4) 

Z jik 



(4) 

^ kij ■ 



(58) 



In these expressions, one has set 



(2) 

Zi = 


1 »w 


(A) 


„3 

'.1 






(3) 

Zij = 


1 a^ 


-.(l+ 7 )^e jmk S% 


<5jfc o-^S^S 
.4 r l . 


"(1 + 7)E^ 

(A) 


(4) 

Zijk = 


~2 di ^ k = 


1G 2 M| 
2 c 4 


r D S 
r 5 


r 5 


D5 t&s&sDs 

r% r% 



GM A D\ 



In the following, let us put 



B. Isotropic geodesic solutions 



T = x" = ct. 



One has now to solve the geodesic equation, written in the following (non afhnc) form 



d 2 x k 
dT 2 



■p/>' 



dx k \ dx a dx p _ 



af} a(i dT I dT dT 



with the isotropic condition 



ffa/3" 



dx a dx p 



dT dT 
Let us first consider the geodesic equation. It develops as 



0. 



■2„fc 



d x 

It 2 



(2) 

r 



_ W* (3)° dx k (3) fe dxl (2)° fai dx k 

oo- r 00 + r 00 dT -2T 0i — + 2T 0idT dT 



(i) dx l dx k (2) fc fai dx 3 (4) fc dx i dx j (3)° da ,i da j dx fe 

+ ai ~dT~dT~ lj ~dT~df~ ij ~dT~dT + ij dT dT If 



The solutions can be written under the form 
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(59) 



(2-4)" 



x k = x k + N k T + x 
where x k is the position at t = (emission). {N k } is any set of three constants related by 

N k N k = 1. 



(60) 



(61) 



(2—4) 

The upper symbol (2 — 4) in x means this part of the solution contains all the (numerically required) cr 2 , c~ 3 

(2) fe 

and c contributions. When useful, the notation x will represent a part of the solution that is only required to 
contain all the c -2 contributions. 

Inserting the connexion (|55|) and (|60p in the geodesic equation (|59"|) . one gets 

sk 



d* i2 x 4) 
dT 2 



(2) , .(2) 



(l +1 )N k J2l3 l Ae Z(A)i 

(A) 



2N l 



(3) (3) . (3) 

Z ik -Z ki +N 3 N h Z i 



(2) 



(2 



.(-1) 



(4) 



S 2 J df 



,@ d 



(2)" 



(2) J 



+2V 00 Z fc + 2N i N j Z ijk - N l WZ kl] + 4N k Z ,•-— + 4iV* Z 



dT 



r/T 



( 2 ) ■ dV 



dT 



The expressions of D^ to consider in this equation read, setting D l Ae = x e — x l Ae and /3 Ae = v Ae /c (values at t — 0) 



Dji (*,**) = ^ e + (iV / -A e )T + 



(2) 



the 



terms 



D\ (t, cc*) = D l Ae + N l T in the c~ 3 and c~ 4 terms 

since, for a hierarchized system, all the /3 terms are at least of order 0(/z) (see section fVI A[) . Now, let us define 

p A = (T + T Ae f-{ 
where (with K Ae > 0) 

T Ae = N l D Ae ~f3 Ae D Ae + 2N l D Ae N m f3? 



K\ e 



Ae " ^/le^ie T iiv ^ e iv p Ae 

1 = (1 + 2N l j3 l Ae ) D% e D% e - T 2 e 



(62) 



= (1 + 27V'/^ e ) £>£*>& " (1 + 4N m ffi e ) (N l D l Ae 
Considering only numerically relevant terms, it turns that 



2ffi e DZN l D A 



r A = {l- N l p l Ae ) 



p\ + 2 (D l Ae + N l T) 



(2)' 



In the c 2 terms, one develops D A /r A and gets 



-^ = (l + 3N m p A l e ) Ae 



D Ae + (N* - (3 Ae ) T 



A 



(2)' 
X 



3 (D> + N l T) { x ] ^ 



N'T 



'Ae 



Pa 



This way, the geodesic's equation now reads, discarding non-relevant terms 
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(63) 



,(2-4)" 



dT 2 



^- - (1 + 7) (S lk 2N*N k ) E^ C 1 + 3 ^ m ^e) 

(A) C 

pa 



^e + (^ - &J ^ 



pi 



_ ( i + 7) ^^df 



-(1 + 7) 


f 1 


G 2 M| 
c 4 


Ps 



Pi\ 



2£ A — 3-^ - 2N*f3 A ^- + 2N k N ] l3 A ^ 

Pa 



4N l e imk ^r + 6e km i- 
Ps 



Pa 
N l D l Q D 



Ps 



n k d 1 
SN% s s 



(64) 



t-iml ' 



Ps 



6N 3 N k , 



N l D l Q D 



i 
s^s 



jml 



GM, 



2N k — ^ - 2^ ^ 



Ps 



Ps 



-AN 



(2) 

k D l s dx 



Ps 



-4- 



■ ■ (2)' 



Ps 



dT 



, (2) 

D%,dx 

2—^N l 



Ps 



dT 



-2 (Sik - 2N*N k ) 



GM S 



{ i-3 { D 7e+ N- T fr£k±Ei 

Ps Ps 



Pi 



Let us point out the development (|63|) is justified only if one has 
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T 



(2)' 



(DL + N l T) x 



Pa 



<3 



/; ; \ (2) 

(D^ e + 7v'r)V 

Pa 



on the whole photon's orbit. In fact, since only the solar terms are required for the higest order contributing terms, 
we consider this condition with respect to the Sun only. Besides, considering applications aimed in this paper, for 
which the contribution of the c -4 terms is expected to be close to the limit of detectability, it is sufficient to require 



{D l Se + N l T) [ x ] 



< gP| X / 



(65) 



where / ~ 1/10, or even 1/100, to keep a safety margin. 

One first obtains the solution for the c -2 and c -3 contributions. Inserting the c~ 2 solution in the rhs of (|6~4"|) . one 
gets the c -4 contribution. Since one considers isotropic geodesies only, the solution has also to satisfy the isotropic 
condition, that writes, using (|57j) and (|60|) (and writting 7 = 1 in the c~ 4 terms) 



(l + ,) { l\o + 2N^\ t+ 2N^ + 2 { l\ N^ 



dT 



dT 



(2)' (2)' 

dx dx 
~dT~~dT~ 



N l N j ^. 



Finally, after some tedious calculations, one gets the solution 

t^GMa 



x k {T) 



■1 , 



G 2 Ml 1 
1-2 — s 



K 2 

^Se 



Wr+(H- 7 )^ 

J (A) 



(2) K (2f 

X A (T)-X A (T = 0) 



(1+7)E 



GM A 



(A) 



,k 



(3) (3)' 

Y A (T)-Y A (T = Q) 



+ (1+7)^; 



r (3) fcm 



3 U S 



(3) 



km 



Z (T) - Z (T = 0) 



G 2 M 2 



(4 )« (4) fc 

X (T)-X (T = 0) 



(66) 
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where 


one has se1 










(2) k 

x A 


= D Ae-N*N l D l AepA NHn(pA + T + TAe) 

K Ae 




(3) fe 

y a 


(S km - N k N™) ((3 l Ae D% e + D A J% e ) N l Qk ,„,,„,, 
= ^ 2 p A + /3 Ae \n(p A + T + T Ae ) 

^Ae 




(3) fcm 
Z 


n D k e -N k N l D l A . / m 1\ 
-2 ^ ATe imn ^ |/s D fl6p J 






\^kmi 


\ Ni ( i i) , 4-»kr+r s ; 






(4) fe 

X 


D k Se - N k N l D l Se 
K 2 

^Se 


4 In 77, 4 - - — arctan — h - — jp - 

PS PSe+Tse \ 4/ K Se K Se 4 p| 


4pSe 4p Se p S 
PS ^ e 




N k 


4:K Se ^ P. 
PS 


5+T + T Se / ex T + T Se 
4 arctan 

PSe + Tse V A J K Se 


T + Ts pSe T + Ts , e 

+ . J<Se 2 + 4 

4 Pg ASe PS 





and with p\ e = p\ (T = 0) = Tj e + K Ae = (l + 2N l /3 l Ae ) D^ e D^ e . The last expression is given with the tracking 
coefficient e (=1 in applications) in order to isolate terms coming from the c -4 part of the metric in the c -4 part of 
the geodesic solution. 

The solution (|66[) has been written in such a way that : 



• all the velocity and spin source terms have been brought together in the c 3 contribution (apart from the j3 
terms entering the definitions of T Ae and K Ae in (1G2[) ) : 

• terms parallel to N and orthogonal to N (i.e. 5 km — N k N m , D Ae — N k N l D l Ae and e kmi N % ) appear explicitely. 
Let us also point out that : 

• thanks to the hierarchy, it is sufficient to put 

T Ae = N l D' Ae (67) 

K Ae = D% e D% e - T\ e = D% e DZ - (N l D l Ae ) 2 

(3) fe (3) fcm (4) fe 

in Y A , Z and X , instead of the full definitions (1521 ; 

(3) fc (3) fem 

• the /3 In (...) term in Y A and tkmiS™ D l Se term entering Z S™ are the only terms that are neither parallel nor 
orthogonal to N. 

For geometrical interpretations, let us remark that : 

• in some sense (among others, since K Ae is defined from the initial position, and not from the (virtual) position 
at t = — oo, see the discussion in [34|), K Ae can be interpreted (from (|67p ) as some impact parameter, or minimal 
approach distance, with respect to the body A; 

• to the same extent, the quantity (D Ae — N k N l D l Ae ) /K Ae can be interpreted as a vector of modulous l/K Ae 
that is pointing toward the closest approach point (with respect to the body A) ; 

• to the same extent, the quantity T Ae can be interpreted as the opposite of time needed to reach the point of 
closest approach to the body A from the initial position x k = x k (T = 0) . 

C. Discussion 

In the inner Solar system, the equation (|66[) gives the time-dependent position for a photon starting from the 
position x k at t = 0. The millimetric level expected in forthcoming experiments has required the post-post-newtonian 
level of approximation. The trajectory is characterized by two arbitrary parameters, giving the initial "direction" in 
some sense. These parameters are, for instance, two any components of the "vector" N k , the third being given by the 
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normalization condition (|61|) . (In relevant applications, the requirement the photon has to reach a given target having 
a known orbit, for instance, fixes these two remaining degrees of freedom and the coordinate time of reception.) 
So far, such post-post-newtonian works have been completed in some one gravitationnal source cases : 

• spherical source in GR [35|; 

• non-spherical source in GR [311 ] ; 

• spherical source in a class of parametrized metric, including GR and, to some extent, ST theories [34j |. 

The current work includes in a coherent way the different gravity sources of the solar system at the required numerical 
level. This includes kinetic source terms effects (velocities of the sources, rotation of the Sun), leading to non-zero 
time-space components of the metric (<?oj); and to the related c -3 contributions in the solution (|66[) . 

However, it has been stressed that (|65[) has to be satisfied on the whole trajectory to ensure the validity of the 
proposed geodesic solution. One has pointed K$e can be crudely interpreted as the closest approach distance to the 
Sun ps.min on the photon's orbit, and Tg e as the opposite of time needed to reach the point of closest approach to 
the Sun from the initial position. The condition (|65p leads to 



2GM S 



I (rp , rp Si T + T Se + PS 

PS - PSe + (T + T S e) In 



TSe + PSe 



< \ps x /■ 



For a given orbit (thence a given ps.min), and if ps e ^ PSm, the most severe condition is reached at the closest 
approach, leading to 

/ 2GM S V 1 2/ 
PSe < Ps, m in x n x T- ( 68 ) 

XPS.rninC J 5 

Since ps.min cannot be smaller than the solar radius, it turns the validity of the proposed geodesic solution is ensured 
as soon as 

the last estimation being obtained taking 2//3 ~ 1/100. This ensures the validity of the proposed solution for inner 
solar system links. 

While it could seem natural to claim that any position on the photon's orbit can be regarded as the initial one, a 
glance at the formula (|6"6"|) makes obvious this cannot be the case. Indeed, at great "distance" (i.e. T — > ±00), one 
has 



(4) 

dX D k Se - N k N l D l Se 

- st — - 



dT J K% 



Se 



e\ 1 1 Apse 



V K Se 2 K% 



Se 



(terms that — > when \T\ — > 00) (69) 



where st is the sign of T. On the other hand, it is clear on physical grounds that dx k /dT should go to values that 

don't depend on the chosen initial point on the orbit when T — > ±00, these limits corresponding to "initial" and 

" final" directions of the photon (the difference of these limits giving the so-called " deviation" of the photon induced 

by the gravitational field). It is clear from (|69|) that dx k /dT computed from ([66]) includes a c -4 term that goes to 

infinity for an initial position chosen (on any given orbit) arbitrarily far from the Sun. This is obviously physically 

unacceptable (and incompatible with the fact that c -4 terms should be corrective terms). This term corresponds 

to one of the so-called "enhanced terms" in [34| (also present in [35[). The validity condition (p5|) . rewritten (|68|) . 

prevents to fall in this case. The tracking coefficient e shows the involved term comes from coupling contributions in 

(2)* 
the geodesic equation (|64[) (i.e from (|64j) -r.h.s. terms depending explicitcly on x ), and not from the fourth order 

(4) . . . . I — 1 

terms g a a in the metric, in accordance with |34| . 

Since the validity condition is ensured only for initial conditions at distances not exceding some astronomical units, 
this shows one should be careful in any extention of such analytical post-post-newtonian approaches to both the 
deviation of light and light tranfer problems at the external solar system scale, at least if photons grazing the Sun are 
considered (remark it is not the case for the GAIA mission). 
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VIII. CONCLUSION 



In this paper, one has get the explicit geodesic equation (|66|) relevant for propagation of light at the millimctric 
level in the inner Solar system, and for related problems. These expressions should be useful for experiments like 
[16l - [l9| . One has also explicitly given the condition of validity of this expression, and pointed out how this validity 
condition is related to the method used to solve the geodesic equation at the post-post-newtonian level. 

This work can be extended to missions involving links at the whole Solar system scale. The planetary c~ 4 metric 
contributing terms are expected to be negligible in this case too, since these terms generate effects of order of 10~ 16 s 
at best, i.e. well under the mm level, even with a large safety margin. Hence, the metric (j50152[) is relevant in this 
case too. However, giant planets multipolar terms have to be taken into account in (|53ti56[) . at least at the c -2 level, 
as shown by previous studies (see [36[ for instance). Besides, the validity condition of the method has to be carefully 
checked in this case, especially when one deals with photons grazing the Sun before/after travelling the external Solar 
system as argued in section [VII CI 
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Appendix A: Effects on the metric of a non-rigid rotation 

Recent results coming from Solar seismology suggest that the Solar core rotates faster than the external layers 
28 4301] . In the following toy model we consider the mass density as spherical (p(x) = p(r)). As usual in Solar models 
26| , we make the distinction between three main regions : the core region (r G [0, Rq ~ 0, 2Rg], where Rq is the core 
radius), the radiative region (r G [Rc,Rr ~ 0, 7Rs], where Rr is the tachoclyne radius) and the convection region 
(r G [Rr, Rs]i where R$ is the Sun radius), tt being the angular velocity, we modelize the differential rotation as 
follows 

n(x) = Q(r)k z , (Al) 



with 



.Re J \ Rr. — Re ) \ Rd — Rr 

where 



n(r) = n o (0)] ( -}-) + n R (9) n ( J R ° )+n D (o)n( J R * ) , (A2) 



n(a;) = lVie[0;l] 

= everywhere else, 



and where we modelize the differential rotation by a simple model in accordance with the usual model 27 1 

n B (9) = n B (l + e B cos 2 8), (A3) 

where fl B and e B are constants, B being any of the three previous regions (C, R or D). 

Solar seismology suggests that the radiative region rotates as a solid - meaning £Lr{6) — Qr (ie. er = 0). However, 
since we are interested in testing solar seismology results, let us relax this assumption. Let us write 

W' s (t,x) = W^(t,x) + Wk(t,x) + Wh(t,x), (A4) 

with 

W B (t,x) = G f d 3 X (il B (R) x X) l ^ R \ . ( A5) 

Jb ^ ' \x — X\ 

Wg being the Sun spin part of W% in ([5T|l . Then the solution writes 

'1 , e B \ Mi , e B M? 



W t s (t,x)=4nGj2{^B 

B 

where 



x x 



t 2 c_B m 4 . o 



3 15/ r 3 35 r 5 



(5 cos z 6 - 1) 



(A6) 



M% = Mjv(0, R c ), M% - Afjv(i?c, Rr), M% = M^^, i? s ), 

with 

f y 



M N (X,Y)= R N+z P {R)dR. (A7) 

ix 

First, note that a faster rotation of a rigid core will modify the value of the total angular momentum only. However, 
this value could be affected by a differential rotation as well. 

But, one also may have to consider a term like the last term of the r.h.s. of (|A6[) in (|55p . in order to measure 
possible weak effects due to differential rotations of the different stages of the Sun - then, giving a characteristic way 
to put constraints on such differential rotations, which will be independent of Solar seismology and neutrino detection 
results. In what follows, we will refer to this term as the diffcrcntial-Lense-Thirring term. 

As an illustration, let us consider the following toy model. Assume (1) the density decreases linearly with the 
distance to the center of the Sun (2) the differential rotation is independent of the distance from the center. Now 
consider that the photons - used for the time transfer - graze the Sun (ie. b = aRs , where b is the impact parameter 
and a(> 1) a parameter ideally close to 1). Then the effect of the diffcrcntial-Lense-Thirring term is about a~ 2 e B /ll 
times the usual Lense-Thirring effect. 

However, realistic models of the Sun and its rotation are expected to substantially decrease this effect. A specific 
study, that considers different realistic models of the Sun, should be done in order to clarify this issue. 



